Statistical properties of currents flowing through tunnel junctions 
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This paper presents an overview of the statistical properties arising from the broadness of the 
distribution of tunnel currents in metal-insulator-metal junctions. Experimental current inhomo- 
geneities can be modelled by a lognormal distribution and the size dependence of the tunnel current 
is modified at small sizes by the effect of broad distributions. 
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I. INTRODUCTION: MIM JUNCTIONS AND 
BROAD DISTRIBUTIONS 

Metal-Insulator-Metal (MIM) tunnel junctions have 
been introduced into the physics toolbox four decades 
ago . They have given rise to several landmarks in con- 
densed matter physics such as the Josephson effect and 
the Coulomb blockade. Since the discovery of large room 
temperature Tunnel Magneto- Resistance (TMR) g |), 
MIM junctions have been under intense scrutiny again. 
This paper will summarize recent studies of the disor- 
der effects in MIM junctions, which is one of the aspects 
of MIM junctions' physics. This topic is not new |3J ||] 
but it can be revisited thanks to recent experimental and 
theoretical developments. 

Practioners know how difficult it is to achieve repro- 
ducibility of the conductances of MIM junctions, even 
when the junctions are prepared on the same wafer. This 
is becoming a crucial problem with the prospect of appli- 
cations of TMR to Magnetic Random Access Memories 
and magnetic read heads, which require the conductance 
dispersion to be typically less than 10%. This raises the 
question of whether the conductance irreproducibility is 
a purely technical problem or whether there is something 
more fundamental behind it. 

The observed large dispersion of conductances from 
one junction to another one is statistically unusual and 
this provides a clue on the nature of the problem. Con- 
sider, for instance, a 10 x 10 /jm 2 junction with a typical 
interatomic distance of 0.3 nm, so that the cross-section 
contains n ~ 10 9 atoms. According to the central limit 
theorem, relative fluctuations of an ensemble of n com- 
ponents scale as 1/y/n. Thus, fluctuations of 10% at the 
junction scale would correspond to fluctuations of 3000 
at the atomic scale. This suggests that either the dis- 
tribution of tunnel currents is extremely broad, or that 
fluctuations do not average out as in the central limit 
theorem, or both. 



During the last fifteen years, the importance of such 
broad distributions has emerged in several areas of sta- 
tistical physics related mostly to anomalous diffusion 
H @, H • The paradigm of broad distributions is the Levy 
flight, i.e., random walks in which the length I of the 
free flight has a power law distribution P(l) = cxIq /l 1+a 
(I > Iq) with a diverging second moment (0 < a < 2). 
With such distributions, the variance is infinite, the usual 
central limit theorem does not apply and the relative fluc- 
tuations of a sum of n terms do not decrease with the 
number of terms. This reminds of the large fluctuations 
observed even in large junctions. Moreover, the sum of n 
terms displacements tends to be dominated by a few of 
them which reminds of the infamous 'hot spots', i.e., of 
filamentary like structures carrying most of the current. 
At last, with Levy flights, even the law of large numbers 
can fail to apply, i.e., the typical sum of n terms might 
not be proportional to n even at large n (case a < 1). 
This suggests that the tunnel current might not always 
be proportional to the size of the junction. 

There seems to be a connection between broad distri- 
butions and MIM tunnel junctions. To clarify the matter, 
one needs to know experimentally the distribution of tun- 
nel conductances (Section II). Then one can study the 
consequences of the current distribution (Section HI), in 
particular the scale effects (Section IV). 



II. 



EXPERIMENTAL DISTRIBUTION OF 
TUNNEL CURRENTS 
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Conducting Atomic Force Microscopy (C-AFM) can 
map the tunnel current flowing through an oxide barrier 
Q. In this technique, the conducting tip of an atomic 
force microscope is scanned in contact with the surface 
of the barrier, while a bias voltage between the tip and 
the bottom metallic electrode creates a current flowing 
through the barrier. In this way, one records simultane- 
ously the topography and the tunnel current. The actual 
resolution of the C-AFM is difficult to estimate. Corre- 
lations studies show that current structures smaller than 
1 nm 2 are resolved. There must, however, be some sort 
of convolution by the finite tip size, which produces some 
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FIG. 1: (a) Topography (black = nm, white = 0.2 nm) and 
(b) tunnel current (log scale, black = 40 pA, white = 1 nA) 
images for an Al oxide barrier (Al deposited on Co, Ar + O2 
plasma oxidized, A\O x thickness ~ 1 nm, see |l6f| ). 



underestimation of the current inhomogeneities. 

The topography and current images for a typical oxide 
(A10 x ) barrier are shown in fig. |l[ The topography (0a) 
is very smooth (roughness ~ 0.2 nm), whereas the loga- 
rithm of the current image (|l|b) exhibits a continuum of 
current values going from ~ 40 pA to ~ 1 nA. The local 
I— V characteristics obtained with C-AFM are consistent 
with tunnelling with transmission much less than 1, even 
at the highest current points which are, thus, not pin- 
holes. This type of experiments has been reproduced by 
several groups with similar results O, O, 13, 14|, p"5|. 
They are especially useful for barrier optimization, for 
instance with respect to oxidation [l7j or annealing 

© II- 

Several types of tiny barrier inhomogeneities can be re- 
sponsible for the large current inhomogeneities which arc 
observed in both C-AFM and numerical simulations [|l9| . 
First, due to the amorphous nature of most insulator 
barriers, the metal-oxide interfaces can not be perfectly 
smooth and fluctuations of the barrier thickness are un- 
avoidable. What is of interest here is the 'roughness' of 
the barrier thickness which is usually much smaller than 
the topography roughness || [ll|. Thickness roughness 
of less than 0.1 nm can generate the fluctuations vis- 
ible in fig. 0b. There must also exist inhomogeneities 
of the barrier height. Indeed changing a single atom at 
a metal-oxide interface can induce local barrier changes 
larger than 1 eV ^0|. Statistics on the barrier parame- 
ters obtained from many local C-AFM I~V curves seem 
to indicate that barrier height inhomogeneities play a 
more important role than thickness inhomogeneities pi[ . 
At last, metal oxides are known for containing electron 
traps which are clearly evidenced in noise studies pl| ] 
(see also Section IV) and impurity states (lj, |22|. Fur- 
ther studies (C-AFM, physico-chemistry of oxides, elec- 
tronic and structural simulation, ballistic electron emis- 
sion spectroscopy [p3[...) are needed to clarify the origin 
of current inhomogeneities. 

Even without knowing the origin of the current inho- 
mogeneities, important consequences can be drawn from 
the knowledge of the statistical distribution P(i) of cur- 
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FIG. 2: Probability distribution P(i) of the tunnel current of 
fig. ^b (open circles). The lognormal fit of the high current 
tail of these data (/1 = —2.5, a = 0.83) is shown as solid line 
(the small current tail is distorted by spurious noise). Also 
indicated are the distributions P(I n /n) of currents flowing 
through groups of n shuffled pixels (Section |lV| ), which peak 
around the average current /ii for large n. 



rents. The distribution P(i) presents in many cases 
P, E4[ a lognormal shape (fig. @): 



P(t) 
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exp 
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LN(/i )£ r 2 )(i), (1) 



where /i is a scale parameter and a is a shape param- 
eter, hereafter called the 'disorder strength'. Note that 
standard \ 2 fitting procedures are not adequate to fit 
lognormal distributions because most of the current is 
usually carried by the tail of the distribution while \ 2 
fitting weights heavily the peak of the distribution, which 
does not carry here a significant current and is strongly 
affected by spurious noise. To take into account the large 
value tail more properly, one can, e.g., fit a parabola to 
the large lni branch of lnP(z) vs lni (a lognormal distri- 
bution is a parabola in log- log scale). 

The occurrence of lognormal distribution of tunnel cur- 
rents is not a surprise |2fJ. Suppose, indeed, that the 
distribution P(d) of barrier thickness d is Gaussian with 
mean /id and standard deviation ad- The current i varies 
typically exponentially with d : 



i e 



-d/\ 



(2) 



where io is a current scale and A is the attenuation length 
of the electronic wave functions in the barrier. By defi- 
nition, the exponential of a Gaussian random variable is 
a lognormal random variable. Thus, the current i has a 
lognormal distribution LN(/i,cr 2 ) with /j = lnio — /id/A 
and a = er^/A p5[ . Importantly, the disorder strength 
<7 is ad /A and not Ud/ [id as one might expect naively. 
Thus, a barrier which appears geometrically smooth 
{o'd/fJ-d <§C 1) might be 'rough' (a — crd/\ > 1) from the 
point of view of current statistics, since typically A <C [id 
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(A ~ 0.05 — 0.1 nm, fig ~ 1 — 2 nm), and generate large 
current inhomogeneities. 

More generally, the tunnel current i depends on the 
barrier parameters pb (thickness, height, voltage, ...) typ- 
ically as: 



i = 9b {Pb) exp [f b (p 6 )] 



(3) 



where gt(pb) and fb(pb) vary less strongly than an ex- 
ponential. If pi, presents small Gaussian fluctuations 
of standard deviation a Pb (a Pb <C fi Pb ) around its av- 
erage value fi Pb , then one has p b = fi Pb + ea Pb where 
e is a Gaussian random variable of order 1 (mean 0, 
standard deviation =1). As a Pb <C fJ, Pb , one can write 
hipb) = fb(v Pb ) + ea Pb fl(fi Pb ) and thus 



i = 9b {Pb) exp [f b (jip b )] exp [ea Pb f' b (fj, Pb )] . (4) 

As 9b(pb) varies slowly, one can neglect its fluctua- 
tions. Thus, i appears as the product of a fixed term, 
9b (pb) exp [fb (p Pb )}, by the exponential of a Gaussian 
random variable, ea Pb f' h (ji Vb ) , which is, by definition, a 
lognormal random variable. 

Thus lognormal distributions of tunnel currents emerge 
as the consequence of small Gaussian fluctuations of the 
tunnelling parameters and are a good starting point to 
investigate the tunnel current statistics. 



III. SIMPLE CONSEQUENCES OF THE 
LOGNORMAL MODEL 

The lognormal distribution of currents (eq. (|IJ)) gives 
rise to peculiar statistical properties. If the disorder 
strength is small [a <C 1), the lognormal is close to a 
Gaussian and the usual statistical behaviours, related to 
narrow distributions, appear. On the contrary, if a is 
on the order of 1 or larger, the lognormal distribution is 
broad and it presents a long tail, just as a Levy flight, 
even if, unlike a Levy flight, it has finite average and 
standard deviation. The broadness of the lognormal also 
appears in the fact that the typical current i l (i.e., most 



probable value), 



can be much smaller than the average current fii 



(5) 



(6) 



indicating a large dispersion. 

Fig. |^ illustrates the differences between narrow and 
broad distributions. Fig. [|a represents random values of 
a narrow distribution, a Gaussian of arbitrary mean fi 
and standard deviation a = 2. The X-coordinate may 
represent the position k in a ID tunnel barrier while the 
Y-coordinate may represent the thickness or height of the 
barrier. All values are of the same order of magnitude, 
/i, to within about a. For the Gaussian, the typical value 
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FIG. 3: Differences between narrow and broad distributions. 
Dots represent random values drawn from the distributions 
described in the text; solid lines represent the partial sums 
of these random values (xlO), see Section [flj|. (a) Gaussian 
(narrow) distribution (b) lognormal (broad) distribution. The 
broken line in (b) gives the typical scaling behaviour (see Sec- 
tion 0). 



and the mean are equal. Fig. ||} represents random val- 
ues of a broad distribution, which is the lognormal aris- 
ing from the exponential of the Gaussian values of fig. ||a. 
Fig. ||b may represent tunnel currents. The appearance 
of the fluctuations is now completely different. The pos- 
sible values cover several orders of magnitude. Neither 
the typical value nor the mean, which differ by a factor of 
400, characterize well the range of possible values. Thus, 
practically, the tunnel current through a disordered bar- 
rier is not well characterized by a single value like its 
mean but rather by the full distribution (two parameters 
for a lognormal). 

The large values observed in fig. ^b correspond to 
the hotspots in tunnelling. These are not pinholes: the 
large currents arise from tiny fluctuations of barrier pa- 
rameters (fig. because the exponential dependence 
acts as a 'fluctuation amplifier'. There is a continuum 
of large current values corresponding to sites that have 
nothing qualitatively special, but rather present small 
quantitative fluctuations of the barrier parameters. With 
the lognormal model, one can estimate the current in- 
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homogeneity pq] by calculating the proportion pa of 
the surface (j>a — P(i)di where i a is a parameter) 
carrying the proportion of the average current (i) 
(pi = f™ iP(i)di/(i)). For small a (a < 0.25), half of 
the total current is carried by roughly half of the sites (no 
hot spots). For large cr, the current proportion carried 
by hot spots becomes more and more important. For in- 
stance, for a — 1 (respectively 2), pi = 50% of the total 
current is carried on average by pa = 16% (respectively 
3%) of the sites with highest current. 

The domination of the current by the few largest trans- 
mission sites is indirectly confirmed by the time fluctu- 
ations of the total current. In certain conditions (small 
enough conditions, low temperatures), the current ex- 
hibits strong telegraph noise indicative of electron trap- 
ping and detrapping on a single trap in the barrier 
pl[ |27], H^] . The large effect of a single electron trap on 
the total current flowing through a junction is a strong 
clue for the dominant role of a few hot spots. 

At last, we comment eq. (|^) for the average current. 
One should first be cautious that, contrary to what is 
frequently assumed j|, ||] , the average [ii is not in general 
what is measured on a single tunnel junctions (see Sec- 
tion IV). For a perfect barrier (cr = 0), we recover the 
current e^ 1 without inhomogeneities, as expected. The 
inhomogeneities generate a correcting 'disorder term', 
e a * . This term is always larger than 1: the barrier in- 
homogeneities always increase the average current, due 
to the non-linear dependence of the current on the fluc- 
tuating parameters. Thus, the average current \ii flowing 
through an inhomogeneous barrier of given average thick- 
ness corresponds to the current i{d c s) flowing through a 
thinner effective homogeneous barrier (d c ff < Pd), as al- 
ready noticed in ||. Yet, in the model based on thickness 
fluctuations, as fi = lnio — aWA, the average current fXi 
still depends exponentially on the average thickness /x^, 
just as the current for a homogeneous barrier. This is 
contrast with the modification of the I — V characteris- 
tics shape by the presence of disorder [f| . 



IV. SCALE EFFECTS 

The broad distribution of currents affects the size de- 
pendence of tunnel junctions giving rise to anomalous 
scaling laws. To understand this intuitively, we have plot- 
ted in figs. ||a and ||b (solid lines) the quantities 



D n /n 



E 

fe=i 



d k /n 



and 



J n/n = ^i k /n 

k=l 



(7) 



(8) 



per unit area flowing through a junction of size n. For 
the Gaussian variable d, D n /n is statistically distributed 
around the mean \x and statistically converges to fi as 
l/-y/n when n increases (central limit theorem). At any 
scale, we have D n oc n. For the lognormal variable i, the 
sum /„ behaves completely differently. At small scales, 
/„ takes small values very different from the mean and, 
as n increases, there is a slow upward trend of /„ towards 
the mean \ii : this is the anomalous scaling we investigate 
here. This upward trend is created by the higher proba- 
bility of larger samples (larger n's) to have a very large 
current peak ik that will significantly draw the sum I n 
towards larger values. This effect does not occur with 
narrow distributions like Gaussians for which the largest 
terms in a statistical sample are not large enough to mod- 
ify the sums significantly. 

These scale effects can be studied theoretically ^(i) . 
The problem reduces to finding the distribution of the 
sum I n of n independent currents ik with the same log- 
normal distribution P(i) = LN(/i, a 2 )(i). For moderately 
broad lognormal distributions, we find that I n is also ap- 
proximately lognormally distributed, as LN(/x n , cr 2 ), and 
has a typical value 



nm (l + C 2 /n) 



-3/2 



which represent, physically, the measured quantities at 
scale n. For instance I n /n is proportional to the current 



(9) 



where \ii is the average current (eq. (g)) and C 2 — e a — 1 
is the coefficient of variation. The typical current is thus 
the product of the usual term n/ii by a correction term 
(l + C 2 jri) 3 ^ 2 . For small junctions (n <C C 2 ), the cor- 
rection term is important. Above a characteristic size 
n c = C 2 related to the disorder strength, the correction 
term tends slowly to 1 and the usual behaviour n/x,- re- 
lated to the law of large numbers is recovered. 

The scaling relation eq. (|^) can be tested experimen- 
tally with the current image of fig. ^. For each n, we 
first sum the currents of groups of n pixels to obtain a 
statistical ensemble of values /„ and then construct the 
histograms presented in fig. ||. The histograms' peaks 
give /* (one must not compute the mean but the typ- 
ical value of I n 's: the mean presents no special scaling 
behaviour). Before grouping them, the pixels have been 
spatially randomized to satisfy the condition of statisti- 
cal independence of the i^'s. The result of this procedure 
for 1^/n is plotted in fig. ||. As predicted, I n /n devi- 
ates strongly at small scales from the constant pn that 
one would expect naively and the deviation is well de- 
scribed by eq. ([)[). The maximum deviation is a factor 
2.8 ~ e 3cr2 /2 f or this good quality junction (cr = 0.83). 
For poorer quality junctions, deviations larger than 10 2 
have been observed |24| . 

If one takes into account the spatial correlations exist- 
ing in the barrier by not randomizing the pixels (black 
squares in fig. [IJ), there is still an anomalous scaling of 
In/n and the convergence towards pn is much slower than 
without correlations. Thus the correlations play a cru- 
cial role in the statistical properties of tunnel junctions 
as they combine with the broadness of the current distri- 
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10° 10 1 io 2 io 3 10 4 10 5 10 6 
n (pixel) 

FIG. 4: Size dependence of the typical tunnel current den- 
sity In/n of the fig. [l] junction. Squares (respectively circles) 
corresponds to unshufned pixels (respectively shuffled). The 
solid line is the theoretical size dependence. The pixel area is 
0.15 nm 2 . However, one can not rigorously convert junction 
sizes from pixels to nm 2 because the effective contact area of 
the C-AFM tip is unknown. 

bution to yield typical currents differing from n\ii even 
for relatively large sizes. 

The size dependence predicted by eq. @ can also be 
tested by measuring the currents of many patterned junc- 
tions of different sizes to obtain the typical current at 
these sizes. This has been done for semi-conducting AlAs 
barriers embedded in GaAs ^9|. Again, the typical cur- 
rent per unit size /* jn is found to increase with the junc- 
tion size, in qualitative agreement with eq. (|^). 

In the tunnel junction community, size dependences 
are frequently analyzed in terms of the product R x A 
of the resistance R by the junction area A and one usu- 
ally checks that fix A does not depend on A, which is 
the normal size dependence. However, anomalous size 
dependences have been reported recently fl30| . For 1 nm 
thick AlO x barriers, a significant increase of R x A is 
found, from 60 ilfim 2 for A = 4 /im 2 to 330 fi/im 2 
for A — 80 /im 2 . As R oc 1/I n and A <x n, one has 
R x A oc nj I n - Eq. ([)]) predicts an increase of I n /n with 
n and thus, one expects intuitively a decrease of R x A 
with A. Therefore, the results of p(| apparently jeop- 
ardize eq. (||). However, the correct theory |2(J predicts 
that both the typical I n /n and the typical n/I n oc Rx A 
increase with the junction size: 

(Rx Af oc (l + CVn)- 1 / 2 , (10) 



which shows how counter-intuitive broad distributions 
can be. 

V. CONCLUSIONS AND OVERVIEW 

Several types of experiments (conductive AFM, noise 
studies, scaling studies) provide a body of evidence that 
the distribution of tunnel currents flowing through MIM 
junctions is broad, even in good junctions. This is a nat- 
ural consequence of the exponential dependence of quan- 
tum tunnelling with the parameters. With such broad 
distributions, the typical value is much smaller than the 
average value so that tunnel currents should not be char- 
acterized by a single number but instead by the full sta- 
tistical distribution. The lognormal distribution is found 
to fit the experimental data in many cases. The shape 
parameter a of lognormal distributions is a convenient 
figure of merit to compare the quality of different junc- 
tions. 

The broad character of the current distribution has 
several implications, some obvious, some less obvious. 
First, the current flows heterogeneously through the junc- 
tion, in a way that the lognormal model can quantify. 
Second, large spatial variations of the current imply large 
time variations, i.e., large noise. Third, the average cur- 

2 

rent varies strongly, typically as e a with the strength 
cr of the disorder. Fourth, the size dependence of the 
typical properties of tunnel junctions (resistance or con- 
ductance) is affected by the disorder, especially at small 
scales. One recovers the usual size dependences at large 
scales but the transition from small scale behaviour to 
large scale behaviour is slow. This transition is further 
slowed down by spatial correlations. Correspondingly, 
the large current inhomogeneities that exist at small scale 
average out slowly when increasing the size of the system. 
Thus, even relatively large junctions exhibit large disper- 
sions of conductances which are the relics of the poorly 
averaged small scale inhomogeneities. 

For applications, it is worth mentioning that the effects 
of disorder increase rapidly when decreasing the junction 
size below a characteristic size related to the disorder 
strength and to the spatial correlations. To achieve bet- 
ter reproducibility, apart from the obvious reduction of 
the barrier disorder, one could also aim at reducing the 
spatial correlations. 

We thank K. Ounadjela who gave the impetus for our stud- 
ies of disorder in tunnel junctions and M. Chshiev, Y. Henry, 
A. Iovan, M. Kelly and D. Stoeffler for discussions. This 
work was partially supported by project 'NanoMEM' IST- 
1999-13471. 
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